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I. INTRODUCTION 

Nonlinear sigma models (NLSM) are of great inter- 
est since they appear as effective models for many quan- 
tum systems in various branches of physics, ranging from 
solid state physics (e.g. quantum Hall effect) to particle 
physics (e.g. theory of light mesons) [THE]. They are 
particularly useful in the description of theories whose 
symmetries are broken below a certain scale. From a 
perturbative point of view, NLSM are considered as fun- 
damental theories only in two dimensions, while in higher 
dimensions they are mainly used as effective theories, 
since the coupling constant has negative mass dimen- 
sion and the theories are not perturbatively renormal- 
izable. However, there is the possibility that in d = 3 
dimensions these models are non-perturbatively renor- 
malizable, i.e. "asymptotically safe" jjj. The concept 
of asymptotic safety has mainly been applied in order 
to provide a non-perturbative renormalization of gravity 
[TU] , It is known that NLSM and gravity exhibit many 
structural similarities, since both are described by non- 
polynomial interactions and have the same power count- 
ing behavior. NLSM can therefore serve as interesting 
laboratory to test and develop techniques for the more 
involved theory of gravity. 

A promising approach to explore the non-perturbative 
quantization of asymptotically safe theories is provided 
by the "functional renormalization group" (FRG) [TOl 
ITT] . It describes the scale dependence of a quantum field 
theory and in particular of its effective action by means 
of an exact functional equation. In order to treat theo- 
ries with nontrivial target spaces like the NLSM, some 
advanced background techniques have been developed 
within the FRG framework, which still require a deeper 
understanding. Nonlinear O(N) models provide a useful 
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testing ground for these conceptual issues, because they 
attracted a lot of attention within statistical field theory 
and, as a result, their critical properties are well-known. 

The FRG approach is a suitable tool to investigate 
the phase structure of physical systems. A second or- 
der phase transition, for instance, is related to a fixed 
point of the renormalization group flow. While we will 
explicitly study the nonlinear 0(N) model, it is expected 
that they possess the same critical properties as their lin- 
ear counterparts. The hypothesis of universality states 
that two short-range theories with the same dimension 
and the same symmetries belong to the same univer- 
sality class. This expectation is strongly supported by 
many explicit computations for d = 3, see [TJ H2H15j and 
references therein. While Monte-Carlo simulations yield 
the critical exponents within an 0(A r )-invariant formula- 
tion, there is further need of non-perturbative analytical 
studies that are manifestly covariant and can confirm the 
critical properties of the nonlinear models. 

Understanding the phase structure of field theories 
with curved target spaces and advancing the techniques 
developed to deal with such theories within the FRG 
are the main motivations for the present work. We will 
study the nonlinear 0(iV)-models in a background field 
expansion [161 117] and compute the effective action in a 
truncation that includes all terms of fourth order in the 
derivatives (see also [18] for some related computation). 
We will derive the flow in arbitrary dimensions and then 
focus on the interesting case d = 3. We will employ a 
manifestly covariant background field method and try to 
improve the formalism by taking care of the specific scal- 
ing of the fluctuation fields. 

Another important motivation of this article is to pro- 
vide explicit results for a direct comparison between the 
FRG approach and the widely used ab-initio lattice ap- 
proach to field theories. An accompanying analysis of 
the renormalization group flow of various nonlinear 0(N) 
models on the lattice by means of the Monte Carlo Renor- 
malization Group (MCRG) method will be presented 
elsewhere [19] . It will be interesting to see how these two 
non-perturbative methods complement each other with 
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respect to the FP structure and critical properties. 

This article is structured in the following way: Sec- 
tion [IT] is devoted to a covariant description of the 0(N) 
models with curved target spaces and the development 



of the background field expansion. In section III we will 
recall the FRG approach and show how to implement 
the background field method. Our truncation of the ef- 
fective action is discussed in section |IV| and the running 
of the remaining couplings is derived. Section |V| presents 
an analysis of the renormalization flow and the critical 
properties of the model. The conclusions are contained 
in section IVH 



II. GEOMETRY OF THE MODEL 

A field of the O(N) non-linear sigma model is a map p : 
£ — » S N ^ X from spacetime £ to the unit-sphere in R*. 
In the present work we will choose for £ the Euclidean 
space M. d . It is convenient to regard the components p a 
as coordinates on the target space S^ -1 . The field space 
of all these maps is denoted by M = {</?: K d -> S' Ar ~ 1 }. 

The target space is the homogeneous space S N ~ 1 = 
0(N)/0(N - 1) with O(N) as its isometry group. It is 
equipped with a unique Riemannian metric h a b that is 
0(./V)-invariant. The microscopic action for the model is 
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S[<p] = k / d d x h ab (<p)d^ a d^ b , 



(1) 



where £ is a coupling constant. Since the fields p> a have 
the status of coordinates, it is natural to regard them as 
dimensionless, such that ( has mass-dimension [(} = d—2. 

The action is invariant under arbitrary reparametriza- 
tions <p — > p'(p) of the fields, provided that the metric 
h a b transforms as a symmetric 2-tensor. In addition the 
model admits the O(N) isometries on the sphere as sym- 
metries. These isometries are generated by vector fields 
Kf(<p) which satisfy a generalized angular momentum 
algebra, 



[Ki, Kj] — —fijtKg , 



(2) 



where fiji are the structure constants of the Lie alge- 
bra of the rotation group. The infinitesimal symmetries 
generated by the Ki are non-linear 



(3) 



From the invariant metric h a b on the sphere one ob- 
tains the unique Levi-Civita connection r a b c and the cor- 



responding Riemann tensor R a 



bed 



h nr h. 



bil 



h„dh 



ad'lbci 



Ricci tensor R a b = (N — 2)h a b and scalar curvature 
R = (N-l)(N-2). 

The Levi-Civita connection on the sphere is used to 
construct 0(iV)-covariant derivatives of the pull-backs of 
tensors on the sphere. For example, given a pull-back of 
a vector on the sphere, its covariant derivative is 



d^ b T b a c v c 



(4) 



The pull-back covariant derivative V will be used exten- 
sively throughout this work. Note that d^p 11 transform as 
vectors, while the coordinates p a are scalars. We also de- 
fine the square of the covariant derivative □ = 5^"V v 
and its Laplacian A = — □. 

In order to construct expansions of invariant function- 
als like ([T]) such that the non-linear symmetries ^ are 
maintained, we will employ a covariant background field 
expansion. We therefore promote the metric h a b(ip) on 
the sphere to a metric h a b(<p) on field space A4 where triv- 
ial spacetime indices have been suppressed for brevity. In 
a similar manner, the Levi-Civita connection, the curva- 
ture tensors and the Laplacian can be promoted to M. 
as well [2U] . It is emphasized that the expansion variable 
ought to possess well-defined transformation properties 
both in the background field p a and in the field </> a . It 
would for instance be a particularly hard task to con- 
struct O(N) covariant functionals in terms of the differ- 
ence 4> a — ip a of two points in field space, as it transforms 
neither like a scalar, nor like a vector under isometries. 

For </> being in a sufficiently small neighbor of tp, there 
exists a unique geodesic in M. connecting tp and (f>. We 
then construct the "exponential map" 



6 a =Exp„r = <£>,0- 



(5) 



Here £ is an implicitly defined vector that belongs to the 
tangent space of M. at ip. Also, £ is a bi-tensor, in the 
sense that it has definite transformation properties under 
both <p a and <p a transformations and for this reason we 
shall write it often as £ a = £, a (<p, <t>). It transforms as a 
vector under the O(N) transformations of <p a and as a 
scalar under those of cf> a . By construction, the norm of 
£ equals the distance between ip and <fi in M. [21] . The 
definite transformation properties make £ a candidate to 
parametrize any expansion of functionals like around 
a background tp. The correspondence is invertible, there- 
fore in (J5J) we introduced the notation <p a = tfi a (<p,{;) to 
be used when <fi a is understood as a function of £ a . One 
always has to keep in mind that the relation between £ a 
and 4> a is nonlinear. 

This background field expansion can easily be applied 
to any functional F [(/>], 



F[4>{<p,t)]=F[<p,€\ 



(6) 



where by abuse of notation we indicated with the same 
symbol the functional when it is a function of the single 
field 0, or of the couple {tp, £}. Even though it may not 
be evident from the way the r.h.s. of Q is written, a 
functional like F[<p] will never be a genuine function of 
two fields, but rather a function of the single combina- 
tion cf> a ((p,{;), and thus is called "single-field" functional. 
There may be, however, functionals of p and £ indepen- 
dently, and those are called "bi-field" functionals. As an 
example, we define the "cutoff- action" 



As k {p 1 a} = 1 - J d d x cK b (^e 



(7) 
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where TZ^ b (ip) is some symmetric 2-tensor operator con- 
structed with ip a . The cutoff action will play an impor- 
tant role in the next section, where we will also specify 
its properties. The functional ^ is still invariant under 
transformations of the background field f a as well as the 
field 4> a ; however, for general TL* h {ip), there is no evident 
way to recast it as a functional of the single field <p a [55] . 

We now construct an expansion of functionals like 
viewed as functions of the pair {f,t;}, in powers of £. 
This can be achieved in a fully covariant way by intro- 
ducing the affine parameter A £ [0, 1] that parametrizes 
the unique geodesic connecting if and (f) [5TJ [53] . Let ip\ 
be this geodesic with tpa = <p and <f\ = <fi. Let us also 
define the tangent vector to the geodesic £a = d<f\/d\ at 
the generic point tf\. We have that ^ is equivalent to 



^1 = T 

dX 

ca £a 

?0 — ? 



/. c(f\ 



(8) 
(9) 



We can rewrite the differential equation satisfied by the 
tangent vectors by introducing the derivative along the 
geodesic Va = £\V a to nn( l * na t = 0. By con- 

struction V \h ab — 0. Another important property is the 
relation to the pull-back connection 



(10) 



The commutator of the pull-back of the covariant deriva- 
tive Va with V M can be computed on the pull-back of a 
general tangent vector v a 



d„.b 



[V A) V> a = i? C( A(^)f «w 



(11) 



and is extensively needed in the covariant expansion. Let 
us now use the covariant derivative along the geodesic to 
perform the expansion of ([6]). Viewing a general func- 
tional F[4>] as limit A — > 1 of F^a] we expand the latter 
in powers of A around A = 0. One can show that 

1 An -i 

F ^ = E ^ *«a=o = E *Ma=o - 



n>0 



ri>0 



(12) 

where we used the fact that F[<f>] is a scalar function of 
<f>. The r.h.s. yields an expansion in powers of £ a of the 
form 



(13) 



n>0 



As an example we give the first few terms of the expan- 
sion of the microscopic action 0, 



S[(j>] = S[<p] +(Jd d x h ab d a f a V^ b 

J d d x (v M rV^ Q + Rabcdd^d^T^ 



2 ' ' "' 

o(e). 



(14) 



III. FUNCTIONAL RG 

We define the scale-dependent average effective action 
of the non-linear O(N) sigma model jTHJ [TT] by the func- 
tional integral 



,-r fc [v>,«] 



x e 



-AS k [ip4-i] 



(15) 



with density of the covariant measure fi[(p] = Det h(ip) 1 ^ 2 . 
Averages are defined through the path integral on the 
r.h.s. of (15), i.e. £ a = (£ a ). Using the relation 
£, a = £, a {fi9) one can obtain 4> a that plays the role of 
full average field of the model. The definition dl5[) differs 



from the usual definitions of effective action only by the 
presence of the functional A5 fc defined in (J7|. This cutoff 
action is chosen in such a way that the average effective 
action Tfc interpolates between the classical action 
at the scale k = oo and the full effective action at k — 0. 
It is chosen to be quadratic in £ a with a ^-dependent ker- 
nel TZK that regularizes the infrared contributions of the 
fluctuations £. This kernel will be called "cutoff kernel" 
from now on. The reference scale k distinguishes between 
infrared energy scales k' < k and ultraviolet ones k' > k 
and the cutoff action leads to a scale dependent effective 
action in the following way: The propagation of infrared 
modes in the path integral (flBl is suppressed, such that 



only the ultraviolet modes are integrated out and we are 
left with a scale-dependent average effective action T k 
for the remaining infrared modes. In order to provide for 
an interpolation between the classical and the effective 
action, the cutoff-kernel has to fulfill the two conditions 



lim fc ^o K„ 



and linn 



KM 



Due to the presence of the cutoff term, is gen- 

uinely a bi-field functional [55J [53] , 



T k [(p,(f>] = r fe 



(16) 



This observation is important in order to understand that 
the only way to construct a single field effective action is 
to set f = 4> or equivalently £ = and consider 



,o]. 



(17) 



The limit k — > of Tk[<f>] coincides with the well known 
effective action introduced by deWitt [25] , 



The very useful feature of the definition (151 is that 



Tk[<f-i £] satisfies a functional renormalization group equa- 
tion QT] 



kd k T k [ip,Z] = -Tr 



-(0,2) 



kd k TZ k [f] 
[<p,£\+K k [<p]), 



(18) 



The functional differential equation ( 18 ) is equivalent to 



the definition (15) once an initial condition La [<£>,£] for 
some big UV-scale A is specified, which accounts for the 
renormalization prescription and the inclusion of counter 
terms. 
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Having an equation like ( 18 1 at our disposal, it is possi- 



ble to investigate properties of the quantum field theory 
under consideration, without having to explicitly com- 
pute the path integral ( |15[ ). In particular we can inves- 
tigate whether or not the theory admits a second order 
phase transition for some value of its coupling constants. 
It is well known that, from a renormalization group per- 
spective, the critical behavior is dictated by a the fixed 
points of the renormalization group flow. 



IV. HIGHER DERIVATIVE MODEL AND BETA 
FUNCTIONS 

It has been stressed in the previous section that the 
flow equation ( 18 ) is non-perturbative in nature such that 



we can explore nontrivial features of the renormalization 
group flow. For this purpose we must find an ansatz for 
the scale-dependent effective action that includes all rele- 
vant operators. We use a covariant expansion of the effec- 
tive action in orders of derivatives and take into account 
all possible operators containing up to four derivatives. 
We furthermore consider the split 



(19) 



where we dropped the overline on the arguments £ and 
4> for notational simplicity. Some comments are in or- 
der: We have already stressed that due to the presence 



of the cutoff action in the definition ( 15 ) the functional 
rfc[y>,£] depends on the two fields {<p,£} separately, and 
not only on the combination 4>(<p, £). In the split ( [l9| wc 
collected the contribution that can actually be written 
as functional of the "full" field <\> into T^l^if, £)], and 
parametrized the rest as r£[^,£], with r^[<p,0] = [2"2"] . 

As an ansatz for r| [4>(ip, £)] , we use the most general 
local action up to fourth order in the derivatives [T71 [55] 



ni<f>] = \ [d d x (ah ab d^ a d^ b + a k h ab nra^ b 



+T abcd {<p)d^ a d> 1 <t> b dv4> c d v <i> d 



(20) 



The isometries of the model are respected only if the 
tensor T abc d is invariant. In the simple case of the O(N)- 
model, there exists a unique (up to normalization) in- 
variant 2-tensor h ab and all higher rank invariant ten- 
sors are constructed from h ab . According to the sym- 



metries T a . 



bed 



T, 



«ab)(cd)) 



that can be deduced trivially 



from ( 20 ) , we see that the most general parametrization 
of T abcd reads 



T, 



bed 



L\ t kh a ( c hd)b 



L2,kh a bh c d ■ 



(21) 



Using this parametrization in (20), wc realize that a total 
of four couplings has been introduced a k , Li.kj ^2,fc}- 
They parametrize the set of operators that we include in 
our truncation and encode the explicit /c-dependence of 

r?. 



Now we need a consistent ansatz for r£ [<£>,£]. It is im- 
portant to dress the 2-point function of the field £ a cor- 
rectly, because it is the second derivative w.r.t. £ which 
determines the flow ((181. Wc choose 



(22) 



+z k n ^ / d d x h ab ee- 



This choice includes a mass term for the fluctuation fields 
as well as a nontrivial wave function renormalization of 

1/2 

these fields £ a — > Z k £ a , which takes into account the 
possibility that the fields ip a and £ a may have a different 
scaling behavior. In a first step beyond the covariant gra- 
dient expansion we assume the the wave function renor- 
malization only depends on the scale k. Later we shall 
see that the wave function renormalization enters the flow 
solely via the anomalous dimension r\ k = —kd k Z k /Z k of 
the fluctuation field. Contrary to Z k the square mass m\ 
enters the flow equation directly. It is the most direct 
manifestation of the fact that is a function of 

the two variables separately. While one could add many 
other covariant operators to T°., we first want to investi- 
gate the effects of this simple structures. 



Now we can plug the ansatz ( 19 ) into the flow equa- 
tion (15). Projecting the r.h.s. of the flow equation onto 



the same operators that appear in the ansatz for T k it 
is possible to determine the non-perturbative beta func- 
tions of the model. In order to proceed we must specify 
the cutoff kernel appearing in Q. We want it to be a 
function of <p solely through the Laplacian and it should 
otherwise be proportional to the metric h ab 



K 



ah I 



= Z k h ab R k [A}. 



(23) 



The cutoff is specified through the non-negative "profile 
function" R k [z] which can be regarded as a momentum- 
dependent mass. The choice (23) ensures that we are 



coarse-graining the theory relative to the modes of the 
covariant Laplacian A. In order to serve as a good cutoff 
the function R k [z] has to be monotonic in the variable z 
and such that R k [z] — for z > k 2 . The wave function 



renormalization of the £ a field has been used in ( 23 ) as an 
overall parametrization. It is convenient to compute the 



scale derivative of (23) already at this stage. We obtain 



kd k K k ab [<p] = Z k h ab (kd k R k [A]-riR k [A}) . (24) 

For the sake of convenience we will compute the beta 
functions of the two sets of couplings a k , L\ yk , £2,*;} 
and {Z k ,m\} in two separate steps. We begin the com- 
putation of the flow of T| [4>((p, £)] by considering the limit 
£ — >• of (fT5| . The result is a flow equation of the form 



kd k ri 



kd k K k [<p] 



-TrGfe (kd k R k [A] - r,R k [A]) , (25) 
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where the dependence on Z k enters via the modified prop- 
agator Gfc which is the inverse of Z^T^' 2 ' [(p, 0] + i?fc[A]. 
It shows how the fluctuations £ drive the flow of the cou- 
plings {Cfc, a k , L\ k , L2.k}- The m odified propagator is 
computed from |S| using ( 20|22 1 and reads 



G fe = (J\[A] + E)- 
P fe [A] = a fc A 2 + CfcA- 



■m 2 I- 



i? fe [A] 



(26) 
(27) 
(28) 



The matrices Z?^, G M and D are endomorphisms in the 
tangent space. In the following we will need the explicit 
form of two of them: 



, ab = 25^{a k R acbd -T abcd )d p tp c d p <p d 



D ab — 

-AT acbd d»ip c d v y d 
Dab = -( k Racbddpp c d p ip d - a k R acbd Dip c D(p d 

+ {0-kRacdeRbfg C 



2-R e (a6) fT e gcd) 



d p (p c d p <p d d a <p f d c 



(29) 



Each term in B pv and D consists of at least two deriva- 
tives of the field ip a . This implies that a Taylor expan- 
sion of ( 26 1 in E is possible, because in our truncation 



ansatz we are interested only in terms up to fourth order 
in derivatives [27]. The tensor C p contains three deriva- 
tives of tp a and thus can be ignored in our truncation. 
Thus the expansion reads 



PZ 



■P-^P-^P- 1 



0(d 6 ) . 



Inserting this expansion into ( 25 ) and using the cyclicity 
of the trace we obtain 



kd k Ti[cp] = iTrA[A]-iTrS/ 2 [A] 



1 



TrE 2 / 3 [A] +0(d 6 



where we defined the functions 

kd k R k [z] - i]R k [z] 



P k [z] 1 



(30) 



(31) 



In ( 30 ) we also used the fact that we can commute £ and 
Pfc[A] in the third term in the expansion, because their 
commutator leads to terms of order 0(d e ). 

The traces appearing in (30 1 arc computed using off- 



diagonal heat kernel methods~{28 . To outline the general 
procedure briefly we consider the traces 



Tr(V Ml ...V^/[A]) 



(32) 



which transform as tensors under isometries. We are in- 
terested in the particular cases < r < 4 and /[A] = 
/; [A] for some I = 1,2,3. Introducing the inverse Laplace 
transform ,C _1 [/](s) of f[z], we rewrite (32) as 



/ d8£- 1 \f\(8)Tt(y lll ...V l . 
Jo 



-sA 



)• (33) 



The expansion of the trace Tr (V Ml . . . V^e - ) is ob- 
tained from the off-diagonal heat kernel expansion (the 
case r = yields the trace of the heat kernel itself). 
This expansion is an asymptotic small-s expansion that 
corresponds, for dimensional and covariance reasons, to 
an expansion in powers of the curvature and covariant 
derivative. Thus we have 



B„ 



TrfV V e ~ sA l = V ^i-Mr.w ± 

n=0 y ' 1 



(34) 



where [r] = r if r is even and [r] = r — 1 if r is odd. The 
coefficients « „ contain a number of powers of the 
derivatives of the field that increases with n, thus only 
a finite number of them is needed to compute the traces 
with 0(d i ) accuracy [28]. The final step in evaluating 
(32) is the s-integration. We define 



in, i 



(4 



i r°° 



(35) 



which we will denote as "Q-functionals" that can be 
rewritten (for positive n) as Mellin transform of f n (z) 



Q, 



l 



dzz^Mz] 



(36) 



(47r) d / 2 r[n] J 

With the help of the heat kernel method we now expand 



the right hand side in the flow equation ( 30 1 and obtain 
kd k i 



■.IfcM = ~tr J d d x{~Q § _ 2>1 [V M , V„] 2 

+ ^f W ( BH ^ + ^ B 'V) 2 ) (37) 



Q, 2 D - Q i+13 B p ^D + Q^D 2 } 



whereby the tensor fields B and D are given in ( 29 ) . The 



beta functions for {(fc, a k ,Li, k , £<2,fc} are defined as their 
log-fc-derivatives and are denoted by f3 a , /3l 1 , Pl 2 }- 
We extract them by evaluating both sides of the flow 
equation at £ = 0. With ( 20 ) the left hand side reads 



+ p Ll (d^d^) 2 + h 2 (d^ a d^ a ) 2 ) , (38) 



and a comparison with the right hand side in ( 37 ) yields 
the beta functions 

PC = Ck( N - 2)% 2 + {N- 2)da k Q, +lt2 

+ L 1:k (N + d)Qa +l 2 - L 2 , k {(N - l)d + 2)Q i+12 
Pa = (N — 2)Q d 2 a k . (39) 

The remaining beta functions Pl 1 and /3l 2 are quite long 
and given in the appendix (Al). It suffices to know that 



G 



Phn Pl 2 are quadratic polynomials of L\ and L 2 and 
that the coefficients of these polynomials depend nontriv- 
ially on both the parameters N and d, and on the other 
two couplings. The combined limit {d = 4, = 0,77 = 
0,m k = 0} of {0L!,0l 2 } is known to yield a universal 
result which can be compared with the results of chiral 
perturbation theory in the case S 3 ~ SU(2) [TTIl2"9"] . We 
checked that we can reproduce the perturbative results 
in this limit. 

Let us finally outline the computation of the flow of Z k 
and mi. The simplest setting to perform this is a vertex 



that depends solely on p 2 . Comparing the power p n with 



expansion of the flow (18) in powers of the field £ a . We 



first notice that for a constant background field ip% the 
ansatz for the effective action ( 19 ) reduces to 



r^c^H— J d d x{( k h ab v^ a v^ b + a k h ab aeae 

4 



a k z k R abcd eaz b aet d 



(40) 



In this particular limit the pull-back connection Q be- 
comes trivial: V M = and □ = d 2 . This observa- 
tion is particularly useful, since we can now easily per- 
form the computations in momentum space. We define 



£°(x) = J d d qe lqx £,q and obtain from @ the 2-point 
function for incoming momentum p M : 

rf 2) [<Pc, 0} P ,-p = Z k (a kP 4 + CkP 2 + ml) . (41) 
We also compute its scale derivative 

fc6> fe lf' 2 V,0]p,- P = Z k ((f3 a - V a k )p 4 (42) 

+ {0c ~ vCk)p 2 + {0m- ~ mnl)) ■ 

On the other hand, the quantity kd k T k °' 2 ^ [ip c , 0] p ._ p can 
be computed from ( 18 ) by applying two functional deriva- 
tives w.r.t. £ a , taking the limit <p a = ip® = const, and 
transforming to momentum space. After these manipu- 
lations, the flow equation (18 1 reduces to 



kd k T k ' 2 \tp c ,0] p - p 



(43) 



-Tr/ 2 (aif< 4 V0] p 



2Z, 



-p,i,-i ■ 



The momentum space 4-point vertex function r^. ' 4 ^ \ip e , 0] 
is obtained from ( 40 ) and is relevant in our computation, 



while the 3-point function vanishes for tp a — const. The 
trace that appears in ( 43 1 consists of an internal trace 



on the tangent space of the model, that involves two of 
the four indices of the 4-vertex, and a momentum space 
integral J d d q/(2ir) d . The final result is a long expression 



0,2,4 of (431 with those of (42) and dividing both 



sides by Z k , we can determine the coefficients 



1 



a ~ r\a k = -{N — 2)a k Qd 2 

O 2 



0C -vCk = ^{N-2)£ k Q ii2 - ({dN -d + 2)L 2 , k 
+ {N + d)L hk - (N - 2)da k ^Q i+12 



p m 2-r)m 2 k = —{N~2)d{d + 2)a k Q i+22 
+ ^{N -2)dQ k Q i+12 . 



(44) 



As anticipated, there is no explicit dependence on Z k , 
because it is a redundant parameter. 

One interesting feature of the method arises at this 
point: Using (39), we can solve the system of equations 
(44) in terms of the two unknown quantities {77, /3 m 2}. 
For a solution to exist, one equation of (44) must be 



redundant and it is a nontrivial check of our computation, 
at this stage, that this actually holds true. The final 
result for the anomalous scaling reads 



V =-{N-2)Q ii2 . (45) 



V. 0(N) PHASE DIAGRAM 

We will now analyze in more detail the structure of 
the /3-functions and the resulting phase diagram. For 
this purpose we focus on three spacetime dimensions as 
a particularly interesting case which has been intensively 
studied [TJ[T11[T5]. Concerning the two dimensional case, 
we just want to mention that our computation repro- 
duces the well-known statement that the theory has no 
nontrivial fixed point, but rather is asymptotically free 

In order to evaluate the Q-functionals and hence the 
explicit running of the couplings in three dimensions, we 
have to choose a specific regulator (p3|) that fulfills the re- 



quirements described in section [III| We decided to choose 
a modified "optimized cutoff" [3Tj : 

Rk[z] = {( k {k 2 -z) + « fe (fc 4 - z 2 )) 0(fc 2 - z) (46) 

Note that the /c-subscript of the couplings will be sup- 
pressed in the following. This choice of regulator enables 
us to calculate the Q-functionals explicitly, 



Or. 



k 2n+2 / (2n + 2 - r) + d t )C 
(4ir) d / 2 r(n) \n(n + l){(k 2 + afc 4 + m 2 ) 1 
2k 2 {2n + 4-r] + dt)a 



n{n + 2){(k 2 + ak 4 + m 2 ) 1 



(47) 



but it renders the system of differential equations rather 
involved, since it has to be solved for the derivatives 
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0.010- 



0.005 



43 0.000 




of the dimensionless coupling a reads in d = 3: 
fin = a + (N - 2) Q 3 2 a 



(48) 



-0.010 

0.00 0.02 



FIG. 1: The flow of the couplings for the truncation with two 
couplings a and £ for TV = 3. The arrows point toward the 
ultraviolet. The removed region lies beyond an unphysical 
singularity which is introduced by the choice of the cutoff. 



dta = kdka — p a and dt( = which also appear on 
the r.h.s. of the flow equation, in order to obtain the run- 
ning of the couplings {£, a, L\, L2}. To study the criti- 
cal behavior we extract the canonical dimensions of the 
couplings and rewrite the equations in terms of dimen- 
sionless couplings £ = fc 2 a = at, L\ = k A ~ d L\ 
and Li = k 4 L*2- The result are the beta functions 
Pa, Pli > @L^' wn i cn are involved rational functions 
and hence not given here explicitly. 

Now we are ready to study the phase diagrams and 
the critical properties that arise from these flow equa- 
tions. We will proceed in a systematic way, by including 
more and more operators in our truncation. The sim- 
plest truncation that only contains the coupling £ was 
already studied by Codello and Percacci [TB] . Their work 
points to the existence of a nontrivial fixed point in di- 
mensions larger than two and hence to the possibility of 
non-perturbative renormalizability. 

We want to add the coupling a and the related 
fourth-order operator. The corresponding renormaliza- 
tion group flow is depicted in Figure [T] (for the case 
N = 3) and confirms the nontrivial fixed point found 
in the simpler truncation. The critical couplings are 
C* = 16(N - 2)/(45tt 2 ) and a* = 0. Note that the ar- 
rows of the flow point into the direction of increasing k, 
i.e. towards the ultraviolet. It is interesting to note that 
the coupling a belongs to an IR-irrelevant operator and 
vanishes at the fixed point. In fact, already the simple 
structure of j3 a that is obtained from (39) reveals that a 



has to vanish for every possible fixed point, as the flow 



Since Q 3 2 is strictly positive for any reasonable regula- 
tor, the only possible fixed point value is a = 0. This 
statement remains true when we include the couplings 
Li,L2- There is also a fixed point for A = 1/a = 0, but 
this is a trivial one whose critical exponents are equal 
the canonical mass dimensions of the operators. It is 
the three-dimensional analogue of the fixed point in four 
dimensions which is discussed in [17] . 

The result a* — agrees with an alternative computa- 
tion of the effective action of the nonlinear 0(iV)-model 
up to fourth order, which is presented in [12j . and in 
which a term oc d 2 (f>d 2 4> is not generated, either. How- 
ever, it is very likely that an extension of the truncation 
to the sixth order in derivatives and an inclusion of oper- 
ators like e.g. Ocj) a D(j> a d fJi (j> b 9 /1 (fib will affect the running 
of a and shift the position of the fixed point. This discus- 
sion will be relevant when we compare our results with 
the renormalization flows obtained by the Monte Carlo 
Renormalization Group (MCRG), that will be studied in 
an upcoming work [T5]. At this point we just mention 
the important fact that both non-perturbative methods 
agree on the structure of the flow diagram, as it is de- 
picted in Figure [T] and hence on the existence of the 
nontrivial fixed point and on the number of relevant di- 
rections at this fixed point. But they differ in the position 
of the fixed point, which has a positive a* in case of the 
MCRG. This is in fact not surprising, since the position of 
the fixed point is not universal and because in the lattice 
calculations the higher order operators affect the flow, 
even though the applied RG procedure keeps track only 
of a truncated operator space, cf. Q15] for more details. 

Since a is not generated in this truncation, the system 
of two couplings effectively reduces to the one-parameter 
truncation that was investigated in detail in [16]. While 
the critical value C* = 16(2V - 2)/(457r 2 ) depends lin- 
early on TV, the critical exponent ^/^| f» is independent 
of N: it is —16/15 for all N. In this sense our computa- 
tion is reminiscent of the one-loop large- N calculations 
[IB] , apart from the small deviation of our critical expo- 
nent from the large- N value —1. It is interesting to note 
that the running of £ in the one-parameter truncation 
of the nonlinear model agrees exactly with the running 
that can be derived by means of the FRG if one regards 
the nonlinear model as the limit of a linear model with 
infinitely steep potential |16j . However, while already a 
simple truncation of the linear model reproduces reason- 
able results for the iV-dependent critical exponents of the 
0(iV)-universality class, this /V-dependence is apparently 
suppressed by taking the limit. 

Let us briefly recall what is known about the critical 
properties of the nonlinear O(N) model in three dimen- 
sions: The second order phase transition of the model is 
described by the nontrivial fixed point and its IR-relevant 
direction. The scaling of the critical coupling is directly 



related to the scaling of the correlation length: 

1 



(49) 



where 0* denotes the eigenvalue of the stability matrix 
evaluated at the fixed point, which corresponds to the rel- 
evant direction. The critical properties of linear and non- 
linear O(N) models are intensively studied, see [I1 H5HI5] . 
and it is generally believed that both theories belong to 
the same universality class. While there are Monte Carlo 
simulations of the nonlinear model which manifestly im- 
plement the nonlincarity of the target space and confirm 
this equivalence [H] , the analytic calculations rely either 
on an explicit breaking of the symmetry and/or on an 
embedding of the target space in a linear space. An im- 
portant motivation for the analytic approach presented 
here is to implement the isometries of the target space in 
every step of the calculation. 

In order to become sensitive to the TV-dependence 
of the critical exponent is, we must include higher or- 
der operators in our manifestly covariant flow equation. 
This agrees with the finding of [16 that the derived /?- 
functions of sigma models with different symmetric target 
spaces coincide in a simple truncation. They differ only 
if one takes higher order operators into account, whose 
number and structure depends strongly on the specific 
type of model. 

Hence we increase the truncation and include the op- 
erator Li(h a i,d^(f> a d l/ (j> b ) 2 . The resulting flow of the cou- 
plings £ and Li of the 0(3) model is depicted in Figure 
[2j whc re we have set the irrelevant coupling a to a* = 0. 
It contains the nontrivial fixed point which was already 
discovered in the leading order truncation and which has 
only one relevant direction. The fixed point exists for all 
N, and while the critical value of L\ is almost indepen- 
dent of N and close to —0.013, the fixed point value 
is an involved expression in N which for N — 3 attains 
the values 0.059. It increases with N such that it ap- 
proaches a linear function with a slope of roughly 0.036 
for large N. Similarly as for the leading order truncation 
the flow diagram has qualitatively the same structure as 
the diagram obtained by the MCRG [19j. The operator 
corresponding to L± is irrelevant in both approaches and 
its critical value L\ is small in comparison to C* • 

Actually there are additional fixed points in the trun- 
cation with coupling £, a and L\, some with negative £* 
and one with quite large values of L\ and (* . These 
could be artifacts of our choice ( 46 ) of the cutoff func- 



tions which may develop singularities for negative cou- 
plings. We could not relate the additional fixed points to 
known critical properties of sigma models and their phys- 
ical relevance remains unclear. We will therefore focus on 
the fixed point depicted in Figures [T] and [2j 

As anticipated, the inclusion of fourth order operators 
renders the exponent v sensitive to the dimension of the 
target space. The TV-dependence of the exponent is de- 
picted in Figure [3j while the numerical values are given in 
Table [TJ The values in the third row denoted by "full sys- 
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-0.020 - 



FIG. 2: The flow of the couplings ( and L\ towards the UV 
for TV = 3. We have set a = a* =0. The fixed point with 
one irrelevant direction is at C,* = 0.059 and L\ = —0.013. 
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FIG. 3: The critical exponent v as function of TV, computed in 
the truncation {(,a,Li}. Depicted are the results of various 
approximations in comparison with average values from the 
literature. 



tern" refer to calculations with the truncation {£, a, 
in which the anomalous scaling 77 of the fluctuation fields 
£ is taken into account. If one sets Z = 1, one obtains 
the values in the second row of Table HI If in addition one 
neglects the fc-derivative of the couplings in kdkR. k on 
the right hand side of the flow equation, that amounts to 
an adiabatic approximation, then we obtain the values in 
the first row of Table [I] At the fixed point the fc-derivative 
of the couplings vanish such that the approximation with 
Z = 1 and the cruder adiabatic approximation yield the 
same fixed point couplings. Since v(N) is a rather in- 
volved and long expression, we tabulated only some se- 
lected values in the rows of Table [T] For a comparison 



9 



Jy 


3 


4 


6 


8 


10 


20 


adiabatic approx. 


0.824 


0.924 


0.969 


0.981 


0.987 


0.995 


with Z = 1 


0.654 


0.756 


0.802 


0.815 


0.820 


0.828 


full system 


0.704 


0.833 


0.895 


0.912 


0.920 


0.931 


literature 


0.710 


0.747 


0.790 


0.830 


0.863 


0.934 



TABLE I: The critical exponent v for various values of N 
and different approximations. The last row contains the best- 
known values from the literature. 



we added the last row which contains the values taken 
from the vast literature about the critical properties of 
the O(N) universality class. For N — 3 and N = 4 
the values are taken from the review pQ, which contains 
the results of many independent computations of which 
we took the non-biased mean values. For N > 5 we 
took the mean values of the results in |151 1321 133j , which 
have been obtained by a high-temperature expansion, a 
strong-coupling expansion and six-loop RG expansion in- 
cluding a Pade-Borel resummation. The corresponding 
values deviate from each other by less than two percent. 

The three truncations of the flow equation presented in 
this work yield a critical exponent v whose iV-dependcnce 
roughly agrees with the results in the literature. The 
values obtained in the adiabatic approximation deviate 
considerably from the references values for small N, but 
show the correct large-TV asymptotic. If one takes the 
running of the couplings in kdkTZ into account, the re- 
sults for small N improve significantly, especially if one 
neglects the wave function renormalization. In this case, 
however, v(N) approaches for large N the value 5/6 in- 
stead of the correct value 1. If in addition one includes 
the wave function renormalization, one obtains a criti- 
cal exponent v(N) which is closer to the reference value 
than the adiabatic result and whose asymptotic behavior 
is better behaved as in the approximation with Z = 1. 
The deviation from the best-known value is maximal for 
N = 5, where it is 14%, and the asymptotic value is 15/16 
instead of 1. This is in fact the value we found in the re- 
duced truncation with just one coupling and agreement 
originates from limAr_ J . 00 (L|/C*) = 0. 

It is not clear to what extent the inclusion of the wave 
function renormalization improves the situation: on one 
hand a running Z improves the asymptotic of v(N) for 
large N and other other hand Z = 1 yields more accu- 
rate results for small N. We included the wave function 
renormalization mainly for conceptual reasons since the 
background and fluctuation fields are treated differently 
in the FRG formalism and hence may possess different 
renormalization properties. This is taken into account by 
admitting a running of Z. 

So far we did not consider the last contribution to the 



functional r|[y>, £] in (22) containing the mass parameter 



and a positive mass parameter. However, the results for 
the critical exponent get worse rather than better and 
are close to the values of the adiabatic approximation. 
This is a surprising finding and certainly requires a better 
understanding. For this purpose, the effects of higher- 
order terms in T\\<p, £] ought to be studied. 

Let us finally include the remaining operator with 
four derivatives Li(h a bd tJi 4> a d ,1 (f> b ) 2 . Although it is of 
the same order as the operators with couplings L\ and 
a, it changes the flow such that there is no nontrivial 
fixed point for the system with couplings {£, a, L\, Li}. 
This statement holds true for all possible modifications 
of the flow, i.e. in the adiabatic approximation, in the 
approximation with Z = 1 and even if we include a 
mass parameter. In fact, already in the cruder trunca- 
tion {£,a,Li} there is no nontrivial fixed point and it 
seems as if the renormalization of the coupling Li is not 
well-balanced. The destabilization of the flow induced by 
the L2-term does not seem to depend on a specific choice 
of regulator. In fact, we used the alternative regulator 
Rk[z] = k d+2 /z and confirmed the existence of nontriv- 
ial fixed points as well as the iV-dependence of v in the 
truncation {£, a, L\}, but the fixed point still disappears 
if one includes Li. 

One may wonder why this term with four derivatives 
destabilizes the renormalization group flow. In the com- 
putation of the full effective action, the renormalization 
of an operator of a given order is always affected by op- 
erators of higher order. These contributions are lost in 
a truncation in which the higher order operators are ne- 
glected. In the present case the b eta f unctions of L\ and 
L2 are quadratic functions, see (All in the appendix, 



mi. We included it in order to examine if terms that 
go beyond the ansatz of a "single-field functional" can 
improve the accuracy of the results. If we consider a 
truncation with couplings {£, a, L\, m 2 } we find the same 
fixed point as before with slightly modified critical values 



and the coefficients of the polynomials must be fine-tuned 
such that both beta functions vanish. We checked that 
the beta function of Li , when evaluated at the fixed point 
for the subsystem consisting of all other couplings, is 
nearly zero. Thus we expect that the inclusion of higher 
order terms will slightly modify the flow in a way that 
one recovers the fixed point and the information about 
the phase transition of the O(N) model, that we already 
detected in the truncation {£, a, Li}. 

However, there could be more subtle explanations why 
the flow of the operator (h a t,d^(f> a d fJ "(j) b ) 2 does not lead to 
a stable fixed point. We only want to mention two possi- 
bilities: There were arguments brought forward recently 
[3~4"] that the regularisation procedure of the functional 
renormalization group given by the introduction of ASk 
may naturally require a corresponding modification of 
the path integral measure, which in turn provides an ad- 
ditional term to the flow equation of the average effective 
action. While this term yields only a renormalization 
of the vacuum energy in theories with linearly realized 
symmetries, it can affect the renormalization of nontriv- 
ial operators in theories whose path integral measure is 
field-dependent. 

The second possibility is that, in order to find a sta- 
ble fixed point for the full system, one has to enlarge the 
truncation as dictated by hidden symmetries involving 
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the background and the fluctuation fields. To this day 
the background field method is the most effective way to 
deal with nontrivial field-space geometries in the frame- 
work of the FRG. Nevertheless, further studies maybe 
needed to understand better which truncations in terms 
of background and fluctuation fields ought to be chosen 
in order to maintain the full reparametrization invariance 
of the theory. An ansatz that is based on the so-called 
Nielsen identities was presented recently in the context 
of gravity [24] and it could be interesting to apply this 
approach to the nonlinear sigma model. 



VI. CONCLUSION 

In this article we applied the functional renormaliza- 
tion technique based on a scale dependent effective ac- 
tion to investigate the renormalization group flow of the 
nonlinear 0(A^)-model. The flow is formulated in a mani- 
festly reparametrization invariant way, so that the results 
do neither depend on any specific choice of coordinates 
on the target sphere, nor on an implicit embedding of the 
nonlinear model into the linear one. Since the symme- 
tries of the theory are realized nonlinearly, we adopted a 
geometric formulation where a background (base-point) 
dependence is introduced in order to maintain the co- 
variance of the model. The background field is used to 
construct a quadratic infrared cutoff term for the fluctua- 
tions, whose purpose is to allow us to effectively integrate 
out the ultraviolet modes while simultaneously respect- 
ing the symmetries of the model. The main achievement 
is the construction of a scale-dependent effective action, 
that is O(iV) invariant for both the transformations of 
the background and the quantum field. The consistency 
of our formalism was underlined by the appearance of 
nontrivial relations between the renormalization flow of 
background and fluctuation operators. 

The model has been studied using a truncation ansatz 
that includes all possible covariant operators up to fourth 
order in the derivatives. The beta functions of the theory 
are identified with the derivatives of the couplings w.r.t. 
the logarithm of the RG scale. They provide important 
informations about the phase diagram of the model and 
we concentrated our discussion in particular on the fixed 
points of the renormalization flow, because they are as- 
sociated to second order phase transitions. The scaling 
of the correlation length is a universal property and is 
directly related to the critical exponent corresponding to 
the relevant direction. Our explicit calculations focused 
on three dimensions, in which the critical exponents of 
the model are well-studied and provide a rich literature 
that can serve as a reference for comparison. 

Our investigation did not immediately point to the ex- 
istence of such universal behavior, because we did not 
find a suitable fixed point for the full fourth order sys- 
tem. However, we discovered that in the restricted sub- 
space of couplings, where one of the higher derivative 



couplings is set to zero (L2 = 0), a fixed point for all 
N emerges. It exhibits one relevant direction, which is 
already present in an one-parameter truncation, and it is 
such that the inclusion of further couplings (L\ and a) 
only adds irrelevant directions. The critical value of the 
coupling a is identically zero in a fourth-order expansion. 
However, the one-parameter truncation is not sensitive to 
the iV-dependence of the critical exponent z/, which re- 
quires higher-order operators. The results we obtain for 
v in the truncation {£, a, L{\ agree qualitatively with the 
pre-existing literature, but show some numerical differ- 
ence that is likely due to the limited truncation ansatz. 
We tested the presence of this fixed point against var- 
ious approximations and choices of the coarse-graining 
scheme, to find that it is a very stable result. 

The renormalization properties of the nonlinear O(N) 
model have also been studied by means of the Monte 
Carlo Renormalization Group and will be presented else- 
where |19[ . The flows that are computed for a fourth- 
order expansion of the model are (apart from the cou- 
pling L2) in good qualitative agreement with the results 
presented here. 

As mentioned, a nontrivial fixed point appears in our 
calculations only if one coupling, L2, is neglected. We 
strongly believe that this is not a pathology of the com- 
putation, but rather it is due to the limited truncation 
considered and the quadratic structure of the computed 
beta functions. However, we also admit the possibility 
that the operator corresponding to the coupling L 2 may 
have a special role in the nature of this phase transition. 
Therefore further investigations are required in order to 
obtain a deeper understanding. 
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Appendix A: Beta functions of L\ and L2 

In this appendix we give the explicit form of the beta 
functions of the couplings L\ and L^. These have been 
omitted from the main text because of their size. To- 



gether with ( 39 1 they complete the system of beta func- 



tions of the couplings of ( 20 ) . Wc computed 
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/3 Ll =[(2iV - 5)ii + 2L 2 - a]Q<j 2 + (2{N + 4) + 4d + d 2 )L 2 + 8L 2 + 4 (rf + 2 )L 2a + d(d + 2)c 



2Li(2(d + 6)L 2 + (d 2 + 3d + 2)a)J Q^ +23 + 2[(d + l)Li + 2i 2 + da]CQ| +lj3 + C 2 Qf ,3 + g<2#-2,i . 



/3 L2 = (d 2 (N - 1) + 2d(7V + 1) + 12)1% + {N + 2d + 6)Lf - 2(d + 2)(2 + d{N - 2))L 2 a 



Qi 



+2,3 



+ 2(d 2 + 2d + 4 + N(d + 2))L 1 L 2 - 2(d + 2)(N - 1 + d)aL 1 + d(d + 2)(JV - 3)a 2 
- 2[L 2 ((JV - 2)d + 2) + L X (JV - 1 + d) - (N - 3)da]CQ^ +h3 + [L x + 2(N - 3)L 2 - {N - 3)a]Q| 2 
+ (TV -3)C 2 Q,, 3 ™ 5^-2,1- ( A1 ) 
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